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Abstract 



We analyze the shot noise in a voltage biased superconducting quantum 
point-contact. Results are presented for the single channel case with arbi- 
trary transmission. In the limit of very low transmission it is found that the 
effective charge, defined from the noise-current ratio, exhibits a step-like be- 
havior as a function of voltage with well defined plateaus at integer values of 
the electronic charge. This multiple charge corresponds to the transmitted 
charge in a Multiple Andreev Reflection (MAR) process. This effect gradually 
disappears for increasing transmission due to interference between different 
MAR processes. 
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In the last few years much attention has been paid to the study of shot noise in meso- 
scopic systems |j],|2|]. These time- dependent current fluctuations are a consequence of the 
discreteness of the charge carriers and their measurement can provide information on corre- 
lations and charge of the individual carriers not available in usual conductance experiments. 
In the case of a current I of uncorrelated carriers of charge q, the shot noise reaches its 
maximum value S = 2ql = Sp i SSOn . This result has been recently used for detecting the 
fractional e/3 charge carriers by measuring the ratio S/2I in the fractional quantum Hall 
regime ||. 

In low-transmitting normal-superconducting N-S structures a doubling of the normal 
Poisson noise-current ratio has been predicted due to Andreev processes where twice 
the electron charge is transmitted. In S-N-S or S-I-S structures the situation is far more 
complex. In this case, it is well established that the main processes contributing to the 
current for subgap bias voltages are multiple Andreev reflections (MAR) 0. At a given 
subgap voltage V the current is mainly carried by MAR processes of order n ~ 2A/eV, 
in which a net charge of ne is transferred. One would then expect an increase of the 
noise-current ratio roughly as 1/V for decreasing bias. This qualitative behavior has been 
recently confirmed experimentally by Dieleman et al. || for a S-I-S tunnel junction. These 
authors give an explanation of the observed shot noise enhancement due to MAR within the 
framework of the semi-classical theory of Ref. [^]. 

In large tunnel junctions a quantitative comparison between theory and experiment is 
prevented by unavoidable uncertainties in the junction geometrical structure. On the other 
hand, there has been in recent years a large progress in the fabrication of superconducting 
point contacts with a reduced number of conducting channels with rather controllable trans- 
mission PJTOd. The subgap structure (SGS) in the I-V characteristics of these systems has 
been quantitatively described with high accuracy using fully quantum mechanical theories 



of transport through a single channel p] , |12j| . One would expect that a similar quantitative 
agreement between theory and experiment could be reached also in the case of noise. 

So far the analysis of current fluctuations in a superconducting single channel contact 



has been restricted to a few special cases. Thus, the excess noise (eV ^> A) for a perfect 



transparent channel has been obtained in Ref. [E3|. The noise has also been analyzed for 



zero voltage and arbitrary transmission in Ref. ||14|| , while the perfect transmission and finite 
voltage case has been addressed in Ref. [ |T5|j . 

The aim of this work is to analyze the shot noise for the whole range of transmissions 
and voltages within the same microscopic model used for the calculation of the current 
in Ref. ||12|| . This analysis would allow to answer the question on whether a well defined 
q = S/2I transmitted effective charge can be associated to the processes giving rise to 
the SGS in the current. We shall explicitly show that this is only possible in the low 
transmission regime where the effective charge tends to integer multiples of the electron 
charge n = 1 + Int[2A/eV]. 

We shall consider the case of a superconducting quantum point-contact (SQPC), i.e. 
a short (L <C £ ) mesoscopic constriction between two superconducting electrodes with 
a constant applied bias voltage V. For the range eV ~ A one can neglect the energy 
dependence of the transmission coefficients and all transport properties can be expressed 
as a superposition of independent channel contributions [TjJ. Thus, we will concentrate in 
analyzing a single channel model which can be described by the following Hamiltonian W2 



H{r) = H L + H R + £ {te^)l\\ uCRu + t* e^^ 2 c^c Lr7 ) , (1) 
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where Hl,r are the BCS Hamiltonians for the left and right uncoupled electrodes, 0(r) = 
0o + 2eVr/h is the time-dependent superconducting phase difference, which enters as a 
phase factor in the hopping terms describing electron transfer between both electrodes. The 
normal transmission coefficient, a, can be varied between and 1 as a function of the hopping 



parameter t (see ref. [12| for details). Within this model the current operator is given by 



Hr) = jE K W/2 cL(r)c fiff (r) - t*e-^)/ 2 cL(r)c L(J (r)) . (2) 



The current-noise spectral density is defined as 



S(u, r) = hj dr' e iWT ' < 8I(r + t')5I(t) + 8I(t)8I(t + r') >= % J dr' e iWT ' K(t, t'), (3) 



where 5I(t) = I{t)— < 7(r) > is the time-dependent fluctuations in the current. 

In the case of a voltage biased superconducting contact both < J(r) > and S(u, r) 
contains all the harmonics of the Josephson frequency u>o = 2eV/h, i.e. < J(r) >= 
J2 n In exp[mo;oT]. These quantities can be expressed in terms of non-equilibrium Green 
functions in a superconducting broken symmetry or Nambu representation Gfj~(t,t') and 
G^ + (t,t') where i,j = L,R defined as 



G+r(r,r') 
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and obey the relation G i • (r, r') 

Then, the mean current and the kernel i^(r, r') in the noise spectral density are given 

by 

e 



< I t) > 



K(t, t') 



Tr 



- 2 {Tr p(r)GjF(r, r')t(r')G R ^(r' , r) + £(r)G^(r, r^V^lV, r)- 
^(^^(r.r'^MG^^r) - f(r)G+i (r, r')tV)^+(r', r)] + (r -> r')} , (5) 

where cr^ is the Pauli matrix, Tr denotes the trace in the Nambu space and t is the hopping 
in this representation 
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In the expression (5) for the kernel we have factorized the two body correlation functions 



V 



following a mean field BCS decoupling scheme ||14|| . The problem of evaluating the Fourier 
components of I and S can be reduced to the calculation of the Fourier components of the 
Keldysh Green functions. An efficient algorithm for this evaluation can be found in Ref. [JTSj . 



We shall restrict our attention to the zero-frequency dc component of the noise S = S(0, r) 
at zero temperature. 

Fig. 1 illustrates the behavior of S as a function of V for different values of the normal 
transmission a. For comparison, the dc component of the current is also shown. As can be 



observed, the more noticeable features in the shot noise are: (i) the presence of a strongly 
pronounced subgap (V < 2A) structure, which persists up to transmissions close to one 
(in the dc current this structure is only pronounced for low transmissions), (ii) In the low 
transparency limit the shot noise subgap structure consists of a series of steps at voltages 
eV n = 2A/n (n integer) as in the case of the dc current, (iii) For higher transmissions there 
is a steep increase in the noise at low voltages, (iv) For perfect transmission the shot noise 
is greatly reduced, (v) In the large voltage limit there is an excess noise with respect to the 
normal case. 

Let us start by analyzing the low transmission regime. In this case it turns out that the 



electronic transport is well described by a sequential tunneling picture |17]] and the current 
can be written as the sum of tunneling rates for each multiple Andreev process times the 
transmitted charge, i.e. Io(V) = e^nY n {V) with T n = (2/h) J du R n {uj). The probability 
of an nth-order Andreev process R n is given by |12 



RJuj) 
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p(uo - neV)p{u) ; u 6 [A, neV - A], (7) 



Lt=l 

where p(w) = |cj|/V^ 2 — A 2 is the dimensionless BCS density of states and p(cu) ~ 
A/V A 2 — uj 2 is the Cooper pair creation amplitude. Expression (7) for R n clearly displays 
the different ingredients in a MAR processes, i.e. it is proportional to the initial and final 
density of states, to the probabitlity of creating n — 1 Cooper pairs and to the probability 
of a quasiparticle crossing n times the interface (ct ra ). In this regime the shot noise adopts 
an appealing form in terms of R n 



(8) 



This expression corresponds to the fluctuations of a random variable (the current) having 
a multinomial distribution. This is a generalization of the simple binomial distribution one 
finds in a N-N contact, which leads to the well known a(l — a) behavior []I],[2|,|3J] . 

This analysis shows that in the limit of vanishing transmission only the lowest order 
process with a non-zero probability contributes to the current at a given bias voltage. As a 



consequence the effective charge shows a step-like behavior described by the formula q/e = 
S/2eI = 1 + Int [2A/eV]. The inset in figure 2(b) shows the effective charge for a = 0.01 
together with the step-like function corresponding to the a — > limit. It is also shown in 
this figure the comparison between the tunnel approximation given by Eq. (8) and the exact 
result for a transmission a = 0.1. Notice that even for this small transmission value there 
are deviations from the simple step-like behavior which are increasingly pronounced when 
the voltage is reduced. These deviations are produced by the contribution of more than one 
MAR process at a given voltage. 

As the transmission increases the sequential tunneling picture breaks down due to the 
interference between different MAR processes contributing to the current at any voltage. The 
charge quantization found in the tunnel regime progressively disappears and is eventually 
washed out when approaching perfect transmission. This is illustrated in Fig. 2(a). 

For intermediate transmissions the analysis of the shot noise becomes rather involved not 
only due to the interference between different processes but also to the increasing contribu- 
tion of fluctuations of higher harmonics of the current. On the other hand, in the perfect 
transmission limit, the analysis is again simplified due to the absence of backscattering. 
Within our theory one obtains a simple expression for the shot noise in this limit 



k 2 
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where R n are the multiple Andreev reflection probabilities given by Rn{u) = l\m=o l r ( u; ~~ 



meV)\ 2 ; uj G [neV, (n + l)eV] and r(u) = (u + i\/ A 2 — uj 2 )/A is the Andreev reflection 
amplitude at an N-S interface. This expression can be shown to be equivalent to the result 



of Averin and Iman in Ref . |T5[ . The great reduction of noise that can be observed in Fig. 1 
for perfect transmission is a consequence of having an Andreev reflection probability equal 
to one inside the gap. 

It is also interesting to analyze the large voltage limit, in which the zero-frequency noise 
behaves as limy^oo S = (4e 2 //i)a(l — a)V + S exc , i.e. the shot-noise of a normal contact with 
transmission a plus and "excess noise" S exc . The excess noise as a function of transmission 
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is shown in Fig. 3. S exc has the same physical origin as the excess current (I exc ), which 
arises from the contribution of the first order Andreev process. One should remark that the 
approach to the asymptotic value is much slower for the shot noise than for the current. We 
obtain that at zero temperature S exc is twice the excess noise of a N-S contact |4j] with the 
same transmission. In particular, this relation yields S exc = 2/5eI exc for the perfect ballistic 
case in agreement with Ref. fl3| . 

A word of caution should be said about the validity of this theory in the limit of extremely 
small bias voltage. In this limit there is another energy scale determined by the inelastic 
relaxation rate 77 (which is a small fraction of A) playing a role in the theory (see Ref. |jl4|j ). 
This finite relaxation rate introduces a cut-off in the MAR processes which determines the 



behavior of current |12j and noise when eV <C 77. The precise behavior at V — > thus 
depends on the actual value of 77. 

In conclusion, we have analyzed theoretically the shot noise in a single channel SQPC 
for arbitrary transmission and bias voltage. We have shown that the shot noise can be 
much larger than the Poisson noise (Sp oisson = 2el) due to the occurrence of multiple 
Andreev reflections in which multiple charge quanta are transferred. In the tunnel regime, 
the effective charge q(V) = S/2I shows a step-like behavior which is a signature of the 
coherent transmission of multiple electronic charges. Our results are consistent with the 
available experimental data of Dieleman et al. ||. These authors used large tunnel junctions 
for which a mean transmission of a = 0.17 was estimated. Although a direct comparison 
with theory is difficult due to uncertainties in the junction structure, the experimental results 
for shot noise and the effective charge are in qualitative agreement with the results of our 
Figs. 1-2 for this transmission range. On the other hand, the predictions presented in 
this letter are amenable to direct experimental test using state of the art techniques for 
fabricating atomic-size contacts PJT0|]. In view of the remarkable agreement between theory 
and experiments for the current- volt age characteristics found in these systems, one would 
expect a similar quantitative agreement when measuring current-current fluctuations. 
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FIGURES 

FIG. 1. (a) Current- voltage and (b) noise- voltage characteristics for different transmissions at 
zero temperature. The values of the transmission are the same in both panels. 

FIG. 2. (a) Effective charge as a function of voltage for different transmissions, (b) Comparison 
between the exact result (full line) and the tunnel approximation given by Eq. (8) (open circles) 
for a = 0.1. The inset shows the exact result (full line) for a = 0.01 and the step-like function 
(dashed line) corresponding to the a — > limit. 



FIG. 3. Excess noise as a function of transmission at zero temperature. 
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